Resonant Pair Tunneling in Double Quantum Dots 
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We present exact results on the non-equilibrium current fluctuations for 2 quantum dots in series 
throughout a crossover from non-Fermi liquid to Fermi liquid behavior described by the 2 impurity 
Kondo model. The result corresponds to resonant tunneling of carriers of charge 2e for a critical 
inter-impurity coupling. At low energy scales, the result can be understood from a Fermi liquid 
approach that we develop and use to also study non-equilibrium transport in an alternative double 
dot realization of the 2 impurity Kondo model under current experimental study. 
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Introduction. Measurements of nonequilibrium shot 
noise in current fluctuations in electronic devices [H has 
became a practical tool to probe strongly correlated sys- 
tems with elementary excitations whose charge, e*, pos- 
sibly differs from the electron charge e, the prominent ex- 
amples being the observation of the Cooper-pair charge 
e* = 2e in normal metal-superconductor junctions 2], or 
fractional charges in quantum Hall samples [3] . Remark- 
ably, many-body physics with unusual emergent excita- 
tions arises even when the interactions occur in a single 
point, e.g., in an impurity in a metal, or in a nanoscale 
quantum dot (QD) connecting few leads. Theoretical 
studies of various quantum impurity problems encoun- 
tered the notion of non-Fermi liquid (NFL) behavior, 
with a fully inelastic scattering of an incoming electron 
into an out-going scattering state which does not include 
any single electron component (3] ■ It is of general interest 
to study shot noise is QD systems showing such elusive 
NFL behavior. 

Nontrivial effective charges emerge even for quantum 
impurity problems showing regular Fermi liquid (FL) be- 
havior, for example in the basic single impurity Kondo 
model, realized by a single magnetic (Kondo) QD cou- 
pled to leads, studies of shot noise [5] lead to a predic- 
tion of a universal fractional charge e* = 5e/3 in the 
low temperature regime which was detected experimen- 
tally _7;], reflecting a combination of single electron and 
two-electron backscattering. The crossover which is typ- 
ically addressed in experiment fS] is very rarely under- 
stood theoretically. In this paper we find that a simple 
and yet unusual "noninteracting-like" picture for trans- 
port of particles with effective charge e* — 2e emerges 
along an entire crossover from NFL to FL behavior oc- 
curring in double QDs in series [^, [13, [HI exhibiting the 
physics of the 2-impurity Kondo model (2IKM). 

The simplest 2IKM consists of two impurity spins 
(5l, Sr), coupled to two channels of conduction electrons 
and interacting with each other through an exchange in- 
teraction K. After the standard "unfolding transforma- 
tion" reducing the two spin- 1/2 channels to four 
chiral Dirac fermions, ipia{x), i — 1,2 — L, R, a =t,ii 
X S {— oo,cx3}, the Hamiltonian becomes H = Hq + Hk 



where Hq = J2],a I dxi/j^^^idxi/^ja and 

Hk - Mi^yi'L) ■ Sl + Mi^y^ ■ Sr + kSl ■ Sr,{i) 

where ct(t) is a vector of Pauli matrices acting in 
spin(channel) space. For this model a NFL quantum 
critical point (QCP) was found at K — Kc ~ Tk 
separating a local singlet FL phase dX K > Kc from a 
Kondo screened FL phase a,t K < Kc. However, more 
realistic models containing inter-channel tunneling. 



Hps = VLRii\il^R + H.c. 



(2) 



[or. Hps = ReVLfl(V'^'''V) — IinVi,ij(V'^T^V') with im- 
plicit sum over spin and channel indices] do not show 
a critical point [3, [lB|- The reason for this is that 
Eq. (m results in a relevant perturbation with dimen- 
sion 1/2 at the QCP [lH, leading to an energy scale 
T* cx Tk\vVlr\'^ + [K - KcY/Tk which is finite even 
at i^T = Kc, below which an effective FL theory takes 
over. Here v is the density of states of the conduction 
electrons. This crossover from NFL to FL behavior is re- 
flected in the conductance of double QDs I,li3,[ii|,[i3, 

and in particular geometries, e.g., the series geometry, 
we were able to calculate it exactly [Tlj. However this 
information is not sufficient to uncover the nature of the 
transport. 

In this paper we study the full counting statistics [ISj, 
^ (PCS) for char ge transfer through a series double QD 
along the full NFL to FL crossover. In general, charge 
is transferred in units of e or 2e. A peculiar situation 
occurs at if ^ Kc, where 2e becomes the basic charge 
unit along the full crossover. This striking behavior is 
not captured in a slave-boson mean field calculation [g^l . 

We also derive a local Fermi liquid Hamiltonian gov- 
erning the physics below T* . Using an enhanced un- 
derstanding of this crossover we go beyond previous 
works [m [23, [23II in determining all coupling constants 
in this effective Hamiltonian and obtain a universal the- 
ory depending only on an energy scale, T* , similar to 
Nozieres FL theory (FLT) for the single impurity prob- 
lem 2J|, and on the new FL boundary condition associ- 
ated with the ratio \K — Kc\/ {v\Vlr\Tk) . This approach 
helps to understand the charge 2e carriers. 
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In the geometry proposed by Zarand et al. 17[ , where 
transport proceeds between two leads connected via one 
QD side coupled to a second QD coupled to another 
lead, exact results on the crossover are not available. 
Nonetheless, we use our FLT to calculate universal non- 
equilibrium transport and noise properties at low ener- 
gies when the NFL critical behavior is destabilized by a 
non-zero K — Kc- Our predictions can be probed exper- 
imentally [ist - 

Full counting statistics. We will obtain the full charge 
transfer distribution in a series double QD tuned to the 
2IKM regime, along the crossover from NFL to FL be- 
havior, using the formulation we developed in Ref. 

In terms of abelian bosonization one can write the orig- 
inal free fermion theory with Hk and Hps — > in 

V2 ' 



terms of 8 chiral Majorana fermions xf^ Xi' 



associated with the real {i = 1) and 



imaginary (i = 2) parts of the charge, spin, flavor and 
spin-flavor fermions {A = c, s, /, X); for a definition of 
these fermions, see Ref. [IjJ. Then the free Hamilto- 
nian is Ho[{x'}] = ^J2^j=i I dxx'jd^Xj, where {x'} = 

{x^,xULx^.Xl,Xl,Xlxi\- The Fermi operator 
gives rise to charge e (and no spin) tunneling from left 
to right, and changes Y — {Np — Np)/2 by 1, Ni being 
the total fermion number in lead i = L, R. 

Turning on Hk, the QCP is obtained a,t K ~ Kc from 
the free case by a change in boundary condition (BC) 
occurring only for the first Majorana fermion, xi(0~) = 
— Xi(0+). For energies <C Tk, the leading terms in the 
Hamiltonian describing deviations K — Kc as well as fi- 
nite Vlr can be written [ll[ in a new basis {x}j where 
Xiix) = x'lix) sgn{x) and Xi = = 2,...,8), as 

Hqcp = Ho[{x}] + SHqcp where [2^] 



SHqcp = i ^ XiXi{0)a- 



Here a is a local Majorana fermion, = 1/2, and 
K-Kc 



Al = ci- 



^, A2 = C2\/TkWVlr\, 



(3) 



(4) 



where ci and C2 are constant factors of order 1. Those 
couplings determine two energy scales Af, Aj, and the 
total crossover scale is A^ = X^ + X^ =T* . The operators 
in SHqcp have scaling dimension 1/2, hence they desta- 
bilize the QCP; below the crossover scale T* = Af + A2 
the system flows to FL fixed points whose nature depend 
on the ratio A1/A2. 

By definition the FCS is obtained from the cumu- 
lant generating function x(m) for the probability dis- 
tribution function P{Q) to transfer Q units of charge 
during the waiting time T (which is sent to infinity), 
x(Ai) = HQe'^^^PiQ)- The cumulants (5"Q) can be 



to a formal equivalence of our non-equilibrium formula- 
tion and that of Schiller and Hershfield \QJ} for a single 
QD tuned to the Toulouse hmit, we can borrow directly 
the results of Gogolin and Komnik for the FCS for that 
model [Ssj; translating between the parameters of the 
two models in the limit T* <C Tk, we obtain 



lnx(A*) 

r 



/OO T 
-lln[l+ ^ A„(.).(e^-"-l) 
-00 n 



(5) 



ri=-2 



Here Aiie 
x 



i^i\xi [npil - rLL) + nuil - np)], ^3(6) = 

4e^fI3-"i(l - '^-R): = An\L^R, np ^ {I + e^/'^)~'^ , 

nL,R = np[e'^eV). The presence of one particle as well as 
two particle transport processes in our model is apparent 
from the dependence of the two terms oc (e^''^ — 1) 
and oc (e^^^A' - 1) in Eq. respectively At K ^ Kc, 
giving Al — A-i — 0, Eq. ([5]) is equivalent to the formula 
for the FCS of spinless noninteracting fermions of charge 
2e transmitted though a resonant level of width ^ T* , 
namely the noninteracting formula [3, |2^ is obtained 
from Eq. (O by the replacement 2/x /i, nj^ p —>■ np{e^ 
eV/2) and adding an overall factor of 2. 

The emergence of two particle resonant tunneling at 
K — Kc follows from Eq. ([3]). In this case SHqcp — 
^A2(V4 + ^z)'* ^^'^ form of a Majorana resonant 
level [27|. This operator changes F by ±1, while a nonin- 
teracting resonant level d with SH ~ X2{ipp+tpp)d+Ii.c. 
changes y[= (Np - Nr)/2] by ±1/2. In both models, 
transport is given by processes of even order in A2 , giving 
AY integer in the resonant level model but AY even- 
integer in the 2IKM. In the limit of large V, in either 
non-interacting resonant level model or 2IKM, one can 
do perturbation theory in A? . The T — conductance, 
in this limit, is given by [nf, [ill G oc X^V^, implying 
a second order process. This follows since the first or- 
der tunneling between a fermionic state [either ipf or ipi 
(i = L, R) in the two models] with energy of order eV , 
and the zero energy level (a or d) does not conserve en- 
ergy. In the small V limit of the 2IKM we can understand 
the charge 2e using the FLT developed below. 

The two-particle processes can be probed by looking si- 
multaneously at the current / = e(S^Q) /T and noise S = 
2e^{S^Q)/T. Eq. © gives the T = V = conductance 

G^jL^got where t = A^/A^ = ^^y^^^.^^/cln^y- 



and go — 2e^/h (or setting /i = 1 as in the rest of the 
paper, go = g^/tt). Using Eq. in Fig. [1] we plot shot 
noise S{V) along the crossover from FL {eV <C T*) to 
NFL {eV ^ T*) regimes for various values of K — Kc (de- 
termined by t). Experiments Q extract effective charges 
by fitting shot noise measurements with the formula 



Sfit = 2e*go t{r)[l ~ t{r)]dr , 



found from (<5"Q) - (-j)"|^ lnx(A*)|^^o. 



In fact, due 



(6) 



where t{V) is extracted from nonlinear conductance mea- 
surements t{V) — t{0) — t, and e* is the effective 
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FIG. 1: Shot noise versus voltage for several values of K — Kc 
(determined by t) aX T = 0. Data is fitted in dashed lines 
using Eq. (O. The inset blows up the FL region, with fits 
of e* = e for {t = 0.15,0.5,0.8) and e* = 2e for t = 0.99 
[K — > Kc, red curve). 



charge. For K ^ Kc {t < 1), Sf^t with e* = e gives 
a good fit for sufficiently small V\ see inset of Fig. [T] 
When \K-Kc\ » Tk\vVlr\ {t < 1), this fit with e* = e 
becomes reasonably good along the full crossover; see 
curve with t = 0.15 in Fig. [TJ For K ^ {t ^ I) the 
fit with e* = e works only for an extremely small range 
eV < ^T*ITk\K - Kc\ ^ 0, and, remarkably, the juU 
curve fits with S'flt with e* = 2e; see curve with t = 0.99. 

Fermi liquid theory. While the above results were ob- 
tained formally from a calculation in terms of refermion- 
ization, it is desirable to understand them from an effec- 
tive interacting theory written in terms of the original 
fermions. The reader may skip the technical details and 
go directly to the FL Hamiltonian, Eq. pT|) . 

The crucial observation is that only the linear com- 
bination (AiXi -|- A2X2)/A of the 8 Majorana fermions 
at the QCP participates in this crossover described by 
Eq. ([3]). It can be shown that the effect of SHqcp is 
to modify the BC for this linear combination by a sim- 
ple sign change. In order to write down the FL fixed 
point Hamiltonian, we define a new basis with mod- 
ified BC, {?7}, where -qi = (XiXi + A2X2)/A sgn(a;), 
m = (-A2X1 + ^iX2)/A, and ?7i = Xi, (« = 3,..., 8). 
Hence we can write the Hamiltonian for the FL fixed 
points as HpL = ^^o + SHpL- We expect the lo- 
cal interaction SHp^ to involve uniquely 771 , which is the 
only field participating in the crossover in Eq. ([3]). The 
only candidate for the leading FL operator with scaling 
dimension 2 is 

SHpL ^ >^FLir]idxm\x=o, (7) 

with XpL oc 1/T*. Comparison of the scattering phase 
shift for Eqs. ^ and (O actually gives Afl = 4/r*. 

Equation ([7]) for the FL interaction may be written in 
terms of fields which are simply related to the original 
fermions. The latter satisfy a FL BC parameterized by 



a unitary S-matrix, 

ipia{0~) = Sij[a]ipja{0^) (repeated indices summed). 

(8) 

We will express Hp^ in terms of single particle scattering 
states 'i/ia incoming from channel (lead) i — 1,2 ~ {L, R) 
with spin a = ±1, 

^ja{x) = 0{x)lpja{x) + 9{-x)Sjj>[a]lljj'a{x), (9) 

satisfying ^'^(0"'") = 5'iQ.(0~). In our left moving con- 
vention, the region a; > (x < 0) corresponds to the 
incoming (outgoing) part of the field. 

To find the S-matrix we should relate the BC of the 
?7's, ?7j(0+) = ??.,(0~), {i = 1,8), to the BC of the Via's, 
Eq. ([5]) . The representation of the ipia 's in terms of the 
?7's (or the x"s or x's) is fairly complicated, however 
quadratic forms of those different fermions are linearly 
related. In particular, consider Jm ='■ Mip : where M 
acts in the channel space. It is straightforward to find 
the coefficients cfj such that Jm = Z]ij=i '^ij ■ x'l x'j '■■ 
Now consider Jm{x) at a; = 0+. Using these linear rela- 
tions together with the BC for the x"s (which depends 
on A1/A2), one can find AI' as function of M and A1/A2, 
such that Jm'{0~) = Jm{0^)- This relation between M 
and M' can be used to find the S-matrix, since Eq. ^ im- 
plies M' = s^Ms. Using this scheme, starting with Hps 
with real VpR we obtain: s[a] — cos(2(5) — iasm{2S)T^ , 

where cos 4(5 = A.A ■ The last equation gives 2S modulo 

Ai +A2 

TT. However, under the transformation Vpp —VpR in 
Eq. (HI, we have 2(5 ^ 27r - 2(5 [5^. Thus 

2(5 = arg(Ai +^2). (10) 

6, which is denoted as the phase shift, changes from 
to tt/2 as function of K, and it takes the value of 
7r/4 at if = Kc- This first exact result for the phase 
shift agrees with the numerical results of Jones [l3| . 
To obtain the S-matrix for complex Vpp, one can 
start with a phase rotation of the tp's, rendering VpR 
real [2^, and in the end, transform back to the orig- 
inal basis leads to a rotation of the S-matrix s 
exp(ir3 arg(VLfl)/2)sexp(-ir3 aig{VLp)/2). 

The FL interaction Eq. ([7|) can be explicitly written 
in terms of the the scattering states ^. To achieve this 
formally we (i) switch from quadratic derivative forms of 
the x's, Oki = XkidxXi + XddxXk, in Eq. to quar- 
tic forms, using iXjdxXj + iXj'dxXj' -■ XjXj' ■■ Xj'Xj 
U 7^ (ii) express : XiXj '■ linearly in terms of quadratic 
forms of the ^''s (these linear relations have the same 
coefficients cf-). The result is H^^ = J dx^Udx"^ (in- 
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dices summed), and [sO] 

cos2 {2S)Oii+ sin2 {2S)022 + sm{4S)Oi2 , 
^^Pi^ - ^* L=o' 

On = f{Jl + J%)-i{JL + JRr, 

O22 = (Jl - Jnf- ■■ f^^ea/j^-^^^ij^e^a^M + H.c. :, 

O12 = t J2 ■ *Lea/3*fl,/3*j7e7A-*j5 : +H-C., (11) 

j=L,R 

where Ji =: Ji ^ *if *i (* = L,R) and Cq/j is the 

antisymmetric tensor. This is the main result of this sec- 
tion. It gives an exphcit form of the interactions between 
the scattering states, related to the original electrons by 
Eq. ([9|) . This interaction is weak in the FL regime allow- 
ing to apply perturbation theory in 1/T*. This universal 
FL Hamiltonian follows by strong restrictions due to a 
large symmetry emerging close to the QCP [l^ and lead- 
ing to the simple form of Hqcp in Eq. ([3]). In practice 
the symmetry at the QCP is reduced by marginal and 
irrelevant operators such as the leading irrelevant opera- 
tor dxXiO' (at the QCP); however they will be associated 
with a small parameter \/Tk, and hence are neglected 
at low energies for T* <C Tk- For finite Vlr, K — Kc, or 
intra-lead potential scattering VlV'lV'-L + VriPI^iPr, ad- 
ditional marginal and irrelevant terms are produced at 
the QCP, part of which were present before. However 
close enough to the QCP and starting with a weak cou- 
pling problem, namely for |i^Vi| <^ 1 {i — L,R), 
and \K — Kc\ ^ Tk, those perturbations can be safely 
ignored. 

The emergence of the basic transport charge 2e for 
K = Kc {6 = tt/A) in the series geometry follows at low 
energies as the only term in Eq. (fTTj) which does not con- 
serve the number of ^'i fermions minus the number of 
^-fl fermions is cx (^'[^ec,/5*[^)(*fl^e^5«'M) -I- H.c. It 
converts a spin singlet pair of ^'l fermions into a spin sin- 
glet pair of '^R fermions (and vise versa). At K = Kc the 
scattering states Vl'L(5'fl) are waves propagating freely to 
the right(left), hence the basic transport mechanism is an 
inelastic 2e hacks cattering. 

To calculate transport properties for the Zarand et al. 
double QD geometry [l7| which is under current exper- 
imental study [2^, it is necessary to calculate the sin- 
gle particle Green's function in the 2IKM. Because the 
electron field cannot be expressed in terms of the Majo- 
rana fermions, we have not been able to calculate this 
throughout the crossover, results being necessarily re- 
stricted to the vicinity of the NFL or FL critical points. 
Here wc consider this system at K slightly different than 
Kc in the FL regimes, T,eV < T*, where our FLT 
can be applied, ignoring particle-hole breaking (Vlr. — 
0). The conductance of this system can be expanded 
as kinglet = 9o[sm''di + (ii{{T/T*f + a{eV/T*f}], 
and Gscrccncd = ffo^os^ 5i - 7i{(r/r*)2 ^ a'{eV/T*f}], 
where is a small phase shift associated with marginal 



potential scattering operators. We assume parity sym- 
metry of the device. Using Eq. (fTTj) , after a lengthy but 
straightforward calculation, we determine universal rela- 
tions: /3i = 71 = 0(1), a = a' = 9/107r^. We also calcu- 
late the shot noise in the FL regime, and define effective 
charges (e*/e) = S/2I for the local singlet FL regime 
[K > Kc), and (e*/e)' = S/2{goV - /) in the Kondo 
screened phase {K < Kc), defined in the limit Si 0. 
Using Eq. ^ we obtain (e*/e) = (e*/e)' = 11/9. Our 
predictions should be contrasted with the measurements 
on single QDs with a = 3/2?:^ ^] and e*/e = 5/3 0. 

We thank J. Malecki and Y. Oreg for very helpful dis- 
cussions. This work was supported by NSERC (ES & 
lA) and ClfAR (lA). 
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